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Abstract 



Conway and Radin's "quaquaversal" tiling of R 3 is known to exhibit statistical 
rotational symmetry in the infinite volume limit. A finite patch, however, cannot be 
perfectly isotropic, and we compute the rates at which the anisotropy scales with size. 
In a sample of volume N, tiles appear in 0(iV 1//6 ) distinct orientations. However, the 
orientations are not uniformly populated. A small (0(iV 1//84 )) set of these orientations 
account for the majority of the tiles. Furthermore, these orientations are not uniformly 
distributed on 50(3). Sample averages of functions on 50(3) seem to approach their 
■ ergodic limits as jV -1 / 336 . Since even macroscopic patches of a quaquaversal tiling 

maintain noticable anisotropy, a hypothetical physical quasicrystal whose structure 
Qh' was similar to the quaquaversal tiling could be identified by anisotropic features of its 

r~| ■ electron diffraction pattern. 



1 Introduction and Results 



This paper concerns statistical properties of Conway and Radin's | p R| quaquaversal (QQ) 
tiling of Euclidean 3-space. A quaquaversal tiling is made by copies of a single wedge-shaped 
tile, appearing in an infinite number of distinct orientations. The set of orientations is a 
dense subgroup of 50(3), and the orientations are in fact uniformly distributed on 50(3): 
Given a smooth function / on 50(3), the average value of / on the orientations of tiles in 
a large patch of M 3 is close to the integral of / over 50(3) (with respect to Haar measure). 
This property, called "statistical rotational symmetry" by Radin | JR1| ] , suggests that a large 



quasicrystal based on the quaquaversal tiling should appear round, and its physical properties 
should be isotropic. We call such a material "asymptotically round" . 
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Figure 1: Part of a QQ tiling 



In this paper we address the question of how large such a quasicrystal would have to be 
to make its roundness manifest. We show that the convergence to rotational symmetry is 
polynomial in the sample size, but with an extremely small exponent. Samples with 10 24 tiles 
have pronounced anisotropy, and even samples with 10 100 tiles have noticable anisotropy. In 
short, the asymptotic roundness of the quaquaversal tiling, while mathematically correct, is 
well beyond the range of any conceivable physics. If a physical material were modeled on the 
quaquaversal tiling, it would exhibit considerably more isotropy than a crystal, while still 
being considerably less isotropic than a liquid or a glass. There is some evidence [|R2|| that 
such asymptotically round quasicrystals may indeed exist in nature. 

The slow growth of complexity is not unique to the QQ tiling, but is a typical feature 
of 3 dimensional substitution tilings. Two dimensional substitution tilings have even slower 
growth, with the number of orientations growing only as the logarithm of the sample size. 
We expect that any 3-dimensional substutution tiling that is qualitatively more isotropic 
than the QQ tiling would have to have a very complicated substitution rule. In other words, 
much of the complexity of the distribution of orientations would have to be built in from the 
start, sacrificing the elegant simplicity of the quaquaversal tiling. 

The orientations of a randomly chosen tile within the QQ tiling is described by a biased 
random walk in a dense subgroup of SO (3). Lubotsky, Phillips and Sarnak [|LPS1| , [LPS2 



considered unbiased random walks on such subgroups in their study of optimal distributions 
of points on S 2 , and some of our methods are related to theirs. Similarly, long time behavior 
of unbiased random walks on infinite cyclic subgroups of 5*0(2) have been studied by Su ||Sj. 
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The QQ tiling is constructed as follows (see ||CR|| for details). We begin with a single 
right triangular prism, \/3 units wide, one unit high, and one unit deep, as in Figure 2a. 
This prism is divided into eight congruent pieces, each similar to the original. The front four 
resulting tiles are shown in Figure 2b, the back four in Figure 2c. We then rescale the entire 
pattern by a linear factor of 2, obtaining a cluster of 8 tiles, each congruent to the original. 




Figure 2a: A QQ tile. 
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We can repeat this process of subdivision and rescaling to obtain a cluster of 64 tiles, 
of 512 tiles, and so on. We embed each cluster of size 8 n within the succeeding cluster of 
size 8 n+1 (there are 8 ways to do this). The union of all the clusters is then an infinite tiling, 
typically covering all of R 3 . 

After a "parent" tile is subdivided, the orientations of the "daughters" are closely related 
to that of the parent. We describe the orientation of a tile by a matrix M G 5*0(3), the 
columns of which are unit vectors pointing along fixed features of the tile. There exist fixed 
matrices g 1: . . . , g 8 G 5*0(3) such that the orientation of the i-ih daughter is Mg i: where the 
orientation of the parent is M. After 2 subdivisions, the orientations of the granddaughters 
are of the form Mgigj. After n subdivisions the orientations of the descendant tiles are equal 
to M times words of length n in the generators {ft}. 

We now see why 2-dimensional tilings have extremely slow growth. 50(2) is abelian, so 
the number of words of length n in a fixed set of generators has only polynomial growth in 
n, while the sample size of a substitution tiling is exponential in the number of subdivisions. 
In other words, the number of orientations is logarithmic in the sample size. 50(3) is 
nonabelian, however, so the number of distinct words of length n grows exponentially with 
n, hence as a power of the sample size. The power, however, need not be large. For the 
quaquaversal tiling we will prove 

Theorem 1.1 Let n > 2. 

If n = 2k is even, then the number of distinct orientations in the n-th subdivision of a 
tile is 24(2 fc - 1) = 24(iV 1 / 6 - 1), where N = 8 n is the total number of tiles. 

If n = 2k + 1 is odd, the number of distinct orientations in the n-th subdivision is 
34 x 2 k - 24 = (34/ v / 2)iV 1 / 6 - 24. 

By Theorem 1.1, a macroscopic sample of, say, 10 24 tiles would appear in over 200,000 
different orientations. The orientations do not, however, appear with equal frequency, and 
they are not evenly distributed in 50(3). In a sample of 8 27 ~ 2.4 x 10 24 tiles, a mere 104 
orientations account for 69% of the tiles, while in a sample of 8 120 ~ 2.3 x 10 108 tiles, 1000 
orientations account for slightly over half the tiles. The asymptotic behavior is given by: 

Theorem 1.2 Let p > 1/28, and let F a {n) be the number of distinct orientations needed to 
account for a fraction a of the 8™ tiles in the n-th subdivision. For any fixed a, 

lim 2- pn FJn) = 0. 

n^oo 

That is, 2-"/ 28 F Q (n) is asymptotically smaller than any exponential in n (i.e. any power 
of the volume), so we speak of F a {n) as growing essentially as 2 n//2S , or as the 84th root of 
the volume. 

Finally, we consider the rate at which the orientations approach a uniform distribution 
on 50(3). Let f(M) be a real- valued function of M e 50(3), and let / ra (M) be the average 
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value of / on the 8™ orientations of the n-th order descendants of a tile of orientation M . 
If / is continuous, then f n (M) must converge to fg (3) fdfJ> for every M, where d\i is Haar 
measure on SO(3). However, this convergence may be quite slow: 



Theorem 1.3 There exists a smooth function f : 5*0(3) — > R and a point M e £0(3) such 
that \f n (M) - J so(3) fdn\ > (0.9938)™ for every n. 



We believe this result to be close to sharp: 



Conjecture 1.4 For any continuous function f : £0(3) — > R and any point M e £0(3), 
there exists a constant C such that \ f„(M) - f so{3) fdfx\ < C2~ n / 112 ss 0(0.9938303)™ for 
every n. 



In other words, it appears that functions do converge to their average values exponen- 
tially in the volume, with the exponents bounded away from 0, but there are exponents about 
as small as 1/336. As in all linear systems of exponential decay, the long-time behavior is 
dominated by the smallest exponents, so we expect the anisotropy of a generic function on 
£0(3) to decay roughly as jV -1 / 336 for sufficiently large N. 



The outline of the paper is as follows. In [Section 2| we examine the set of orientations that 



can arise from successive subdivision of a single tile. This set turns out to be a subgroup, 
denoted G(6,4), of £0(3). The analysis of this section is a special case of some general 
results of RRS1| . We establish a canonical form for an element of 0(6, 4) as a word in two 



generators, and provide an explicit list of the orientations in the n-th subdivision, thereby 
proving Theorem We also develop a measure, called "size" , of the complexity of an 



element of G(6, 4). 



In |Section 3| we consider the populations of the various orientations of the descendants 



of a given tile as a function of n, the number of subdivisions. We show that the "size" 
of the orientation of a randomly chosen tile is described by a biased random walk (on the 
non-negative integers) with a reflecting barrier. Although the maximum size is n/2, the 
distribution of sizes is approximately Gaussian with mean n/28 + 0(1). Since there are 
0(2™/ 28 ) words of size n/28 or less, [l'heorem 1.2| follows, 
i i 

Theorem 1.2 is an asymptotic statement. To obtain results for specific values of n 
we numerically implement the random walk for n up to 200. This establishes the claims 
previously made about samples of size 8 27 and 8 120 . 

In |Section 4| , we consider the rate at which the distribution of orientations approaches 
uniformity on £0(3). There is a linear map £ on L 2 (£0(3)) such that /„ = £ n f. £ is similar 
to the Hecke operators considered in | [LPS1| , [LPS2| j , except that those Hecke operators were 



Hermitian, while £ is not. Nonetheless, the spectrum of £ appears to be entirely real. If 
/„ is to approach a limiting function exponentially, the base of the exponent (for generic /) 
is the largest eigenvalue of £ that is less than 1. The operator £ is the direct sum of finite 
dimensional blocks, each block corresponding to an irreducible representation of £0(3). By 
numerically diagonalizing blocks corresponding to the first 300 representations (blocks with 
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total dimension over 36,000,000), we find a large number of eigenvalues over 0.993, with the 
largest being approximately 0.99381, and conjecture that the eigenvalues are bounded above 
by 2- 1 / 112 « 0.9938303. 



2 The group G(6, 4) 



Let R^ denote a rotation about the x axis by an angle 9. The group generated by R 27r / p and 
Hy n / q is denoted G(p,q). (For a discussion of G(p,q) for general p and g, and for further 
generalizations, see ||RS1| , |RS2| .) The set of orientations in the QQ tiling is precisely (7(6,4), 
generated by S := RJ/ and T := R^ 3 . Indeed, the orientations g 1 ,...,g 8 of the eight 
daughters of a tile in standard position all are words in S and T, and it is not hard to get 
S and T themselves as products of the <fc's. To understand these orientations, therefore, we 
must understand (7(6,4). 



A word a in S and T is said to be in canonical form if 

= S a ° T*° 

with a G {0, 1, 2, 3} and &o £ {0, 3}, or if 

a = WST bl ST b2 ■■■ST bn E, 



<\ ,v" / (2.1) 



(2.2) 

with n > 1, with W of the form ( |2.1| ), with each hi = 2 or 4, and with E = 1 or S. We shall 
soon see that each element of (7(6, 4) can be expressed by a unique canonical form. The size 
of an element of (7(6,4) with canonical form fl2.2| ) is defined to be n. (Elements of the form 
( |2.1D have size zero). 



To understand the group (7(6,4), we represent S and T by explicit 3x3 matrices. 
Relative to the standard basis {ei, e%, e%} of R 3 , the matrices of S and T are 



S 




and T 





(2.3) 



The factors of \/3 are inconvenient to work with. We eliminate them by working in the basis 
{ei, v^3e2, 63}, relative to which the matrices of S and T are: 



S 



Three other useful matrices are 






(1 








T = 





1/2 


-1/2 




x0 


3/2 


1/2 



(2.4) 



rp3 



(I 
-1 
\0 








ST 



±2 








±3/2 -1/2' 
-1/2 T l/2 




(2.5) 
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Theorem 2.6 Every element g G (2(6,4) can be expressed in canonical form, and the canon- 
ical form is unique. If g has size n, then, in the representation (2.4), 2 n g is an integer 
matrix with at least one odd matrix element. There exist exactly 8 group elements of size 
and 16 x 2 n elements of size n > 0. 



Proof: Note that 2ST 2 and 2ST 4 are matrices with integer entries of the form 

/ even odd odd \ 

even odd odd , (2.7) 
\ even even even J 

and that the product of two or more matrices of the form given in (|2.7|) is also of the form 
given in ( |2.7D . If g can be expressed in canonical form ( |2.7| ), then 

2 n g = W{2ST bl ){2ST b2 )---{2ST bn )E 



Inform 2.7) (form |J) ■ ■ ■ (form ^7\)E 



= Inform \T$E. (2.8) 

Since S is a permutation matrix (up to sign), and T 3 is the identity (up to sign), 2 n g is 
an integer matrix with four odd matrix elements and five even matrix elements. If g is 
of canonical form (|2.1| ), then 2 n g = 2°g = g is an integer matrix with three odd matrix 
elements. This proves the second statement of the theorem, and also shows that the identity 
can be put in canonical form in only one way, namely 1 = S°T°. 

To see that every element of G(6, 4) can be put in canonical form, begin with an arbitrary 
word 

a = S ai T bl ■ ■ ■ S a "T bn (2.9) 

in S and T. If any of the a^s is even, or any of the biS is a multiple of 3, we can use the 
relations 

ST 3 = T 3 ^ 1 (2.10) 

and 

S 2 T = T^S 2 (2.11) 

to shorten the word. So we can assume that all the a^'s (except possibly ai) are odd, and 
that none of the b^s (except possibly b\) are multiples of 3. We then use ( |2.1(J| ) to extract 
powers of T 3 and move them leftwards, changing T 1 to T 4 and changing T 5 to T 2 . We 
then use ( p.ll| ) to extract powers of S 2 and move them leftwards, changing S* 3 to S. In the 



process, we may change some T 2 's to T 4, s and vice- versa, but that does not affect the form 
( ^2|) . Finally, if b x is not divisible by 3, we use (|2TTTJ| ) to rewrite S ai T bl as S a °T 3b °ST b 'i with 
b[ = 2 or 4. 

To see that a canonical form is unique, suppose g and g' are distinct canonical forms for 
the same group element. Using the relations ( |2.10| ) and ( |2.11| ), we could then shorten the 
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word g'g 1 to a nontrivial word in canonical form, which would have to represent the identity. 
But we have already seen that there are no nontrivial canonical forms for the identity. 

Since group elements are in 1-1 correspondence with canonical forms, we can count 
group elements of any given size. There are 8 choices for W, 2 choices for each fej, and 2 
choices for E, for a total of 16 x 2 n elements of size n > 0. For size there are only 8 choices, 
4 for a and 2 for b . ■ 

We now have the machinery to understand the orientations of tiles in the n-fold subdi- 
vision of a single quaquaversal tile. In a single subdivision, the daughters are described by 
group elements {g±, . . . , g 8 } with 



Note that T = gjg§, so the semigroup generated by the (?j's is precisely the group 
generated by S and T. Also note that only two of the ^'s, namely g 7 and g 8 , have size 1, 
while the rest have size zero. Since g 2 , g 2 , g^gs and g%gi all have size 1 or less, the orientations 
in the 2nd subdivision have size at most 1 and the orientations in the 2/c-th subdivision have 
size at most k. This shows that the number of possible orientations in the n-th subdivision 
is bounded by a constant multiple of 2 n//2 . 

However, we can be more precise. The following theorem implies Theorem 1.1: 

Theorem 2.13 The orientations appearing in subdivisions of a single tile (of standard ori- 
entation) are as follows. Let P m be the set of words of length m in the two generators ST 2 
and ST 4 . For k > 0, the orientations of the 2k-th subdivision are: 

1. T b P, m S a , with no restriction on a or b, and with m ranging from to k — 2, 

2. T b P k ^iS a , with b = 0, 2, 3, or 5 and a arbitrary, 

3. T b P k _ 1 S a with b = 1 or 4 and a = or 2, 

4. ST 2 P k ^S a , with a = or 2, and 

5. T 3 (5T 4 )P fc _ 1( 5 a ; with a = 0or2. 

For k > 0, the orientations in the (2k + l)st subdivision are 



91 = 

92 = 

9s = 

94 = 

95 = 
9e = 

97 = 

98 = 




1 



(2.12) 



= S 2 T 2 
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1. T b P m S a , with no restriction on a orb and m ranging from to k — 1, 

2. T b P k S a , with b = 0, 2, 3 or 5 and a = or 2, 

3. ST 2 P k ^S a , with a = l or 3, and 

4. T z ST 4 P k ^S a , with a = l or 3. 

Proof: It is not hard to see that, in each case, the sets are disjoint and their cardinalities 
add up to the total indicated in Theorem 1.1. To complete the proof we must check i) that 
the claimed orientations for 2 subdivisions are correct, ii) that the claimed orientations for 
2k subdivisions times the generators gi,...,g$ give precisely the claimed orientations for 
2k + 1 subdivisions, and iii) that the claimed orientations for 2k + 1 subdivisions times the 
generators gi, . . . , g$ give precisely the claimed orientations for 2k + 2 subdivisions. We leave 
this straightforward (but quite tedious) computation as an exercise for the highly motivated 
reader. ■ 



3 Population dynamics 



In Section 2] , we counted the number of distinct orientations of tiles in the n-th subdivision 



of a QQ tile, which equals the number of algebraically distinct words of length n in the 
generators g±, . . .g 8 . In this section we consider the distribution of tiles among the various 
orientations. Equivalently, we count repetitions in the words of length n in the g^s. 

We do this by modeling the subdivision process as a random walk on G(Q, 4). If at some 
stage we are at the group element M, then at the next stage we can go to Mg%, Mg 2 , . . . , Mg$, 
each with probability 1/8. The distribution of orientations in the n-th subdivision of a tile 
(with standard orientation) is precisely the n-th iterate of this random walk, starting at the 
origin. The following theorem implies Theorem 1.2 . 



Theorem 3.1 Consider the random walk starting at the origin, and let P(k,n) be the prob- 
ability that the n-th iteration yields a group element of size k or greater. For any p > 1/28, 
lim n ^oo P(np, n) = 0. 



Proof: 

We divide G(6, 4) into sets according to size and the pattern at the end of the canonical 
form. The random walk on G(6,4) then induces a random walk on the space of sets, and 
we show that the bias in this random walk causes the average size to increase by 1/28 each 
turn. The result then follows from the central limit theorem. 

Let R(k, ST 2 ) be the set of group elements of size k whose canonical forms end with 
ST 2 . R(k,ST 4 ), R(k,T 2 S), and R{k,T 4 S) are defined similarly. R(0) is the set of elements 
of size 0. 
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If M e R(k, ST 2 ), then Mg u Mg A and Mg 5 are equal to M, and so are in R(k, ST 2 ). 
Mg 2 is in R(k,T 4 S), Mg 3 is in R(k,T 2 S), and Mg 6 is in R(k,ST 4 ). Mg 7 may be in 
R(k - 1,T 2 S) or R(k - l,T 4 S), depending on whether M ends with ST 2 ST 2 or ST 4 ST 2 . 
However, if Mg-j e — ljT 2 ^), then Mg 8 G — 1,T 4 S), and vice-versa. In either 
case, starting in R(k,ST 2 ), there is a 3/8 probability of staying in R(k,ST 2 ) and a 1/8 
probabilty for each of the destinations R(k, ST 4 ), R(k,T 2 S), R(k,T 4 S), R(k - l,T 2 S) and 
R(k-1,T 4 S). 

This analysis, and the corresponding analysis for R(k,ST 4 ), R(k,T 2 S) and R(k,T 4 S), 
is summarized in the diagram 



R{k-1,T 2 S) 



R(k-1,T 4 S) 




9S or g-j 



gi,94,g5 

n 

R(k, ST 2 




91,94,95,96 

n 

R(k,T 2 S) 



98 or 97 



ST* R(k,ST 4 ) R(k,T 4 S) 

91,94,95,97 91,94,95,96 




R(k+1,ST 2 ) 

97 



R(k+1,ST 4 ) 



This diagram continues to the right indefinitely, but to the left it ends at k — with 

R{1,ST 2 ) 

97 



91,92,93, 
94,95,96 




i?(l, 



Since the long-time behaviors of positively biased random walks are independent of the 
boundary conditions at k — 0, the reflecting barrier at k — does not affect the theorem. To 
establish the theorem, we need only calculate the bias in the random walk away from k — 
and show it equals 1/28. 

We first compute the distribution of ST 2 , ST 4 , T 2 S and T 4 S endings. Summing over 
k, we find the transition matrix for the four classes of group elements is 



(3.2) 





( 3 


2 


1 


3\ 


1 


1 


4 


3 


1 


8 


2 


1 


4 







\2 


1 





4/ 
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The limiting distribution is the eigenvector with eigenvalue 1, namely 



1 



4 
3 



(3.3) 



14 



\3/ 



That is, asymptotically 4/14 of the tiles have orientations whose canonical forms end with 
ST 2 , 4/14 end with ST 4 , 3/14 end with T 2 S and 3/14 end with T 4 S. 



On each iteration, an ST 2 state has a 2/8 chance of decreasing k (and none of increasing) , 
ST 4 states have no chance of either increasing or decreasing k, while T 2 S and T 4 S states 
each have a 2/8 chance of increasing k and no chance of decreasing. Weighing these chances 
by the (4/14,4/14,3/14,3/14) limiting distribution, we see that at each step in the random 
walk there is a 3/28 chance of increasing k and a 2/28 chance of decreasing k, for a net bias 



Theorem 3.1 is an asymptotic result that does not indicate have many iterations are 
needed to reach the asymptotic regime. Indeed, the bias in the random walk is so small that, 
for moderate numbers of iterations, one would expect diffusion and the reflecting barrier at 
k = to be quite significant. To measure this transient behavior we rely on the following 
numerical implementation, in MATLAB, of the biased random walk. Here zz(n) is the 
fraction of orientations of size zero in the n-th generation, st2(k,n) (resp. st4(k,n) , 
t2s(k,n) , t4s(k,n)) is the fraction in R(ST 2 ,k) (resp. R(ST 4 ,k), R(T 2 S,k), R(T 4 S,k)) 
in the n-th generation, and tot(k+l,n) is the fraction of size k in the n-th generation. 

°/ Initialize matrices to zero 

for n=l:Nmax 

zz(n)=0; 

for k=l : (2+Nmax/2) 

st2(k,n)=0; st4(k,n)=0; t2s(k,n)=0; t4s(k,n)=0; tot(k,n)=0; 
end ; end ; 

°/„ Set up distribution for n=l and iterate to n=Nmax. 

zz(l)=6/8; st2(l,l)=l/8; st4(l , l)=l/8 ; tot (1 , l)=6/8 ; tot(2,l)=2/8; 

for n=2:Nmax 

zz(n)=(6*zz(n-l)+2*st2(l,n-l))/8; 

st2(l,n)=(zz(n-l)+3*st2(l,n-l)+2*st4(l,n-l)+2*t4s(l,n-l))/8; 
st4(l,n)=(zz(n-l) + st2(l ,n-l)+4*st4(l ,n-l)+2*t2s(l ,n-l) )/8 ; 
t2s(l,n)=(st2(l,n-l)+st4(l,n-l)+4*t2s(l,n-l)+st2(2,n-l))/8; 
t4s(l,n)=(st2(l,n-l)+st4(l,n-l)+4*t4s(l,n-l)+st2(2,n-l))/8; 
tot(l ,n)=zz(n) ; tot(2,n)=st2(l,n)+st4(l,n)+t2s(l,n)+t4s(l ,n) ; 
for k=2: (l+Nmax/2) 

st2(k,n)=(t2s(k-l,n-l)+t4s(k-l,n-l)+3*st2(k,n-l)+2*st4(k,n-l)+2*t4s(k,n-l))/8; 
st4(k,n)=(t2s(k-l,n-l)+t4s(k-l,n-l) + st2(k,n-l)+4*st4(k,n-l)+2*t2s(k,n-l) )/8; 
t2s(k,n)=(st2(k,n-l)+st4(k,n-l)+4*t2s(k,n-l)+st2(k+l ,n-l))/8; 
t4s(k,n)=(st2(k,n-l)+st4(k,n-l)+4*t4s(k,n-l)+st2(k+l ,n-l))/8; 



of 1/28. ■ 
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Median size = < 



tot(k+l,n)=st2(k,n)+st4(k,n)+t2s(k,n)+t4s(k,n) ; 
end ; end ; 

tot % Output results. 

Running this program with Nmax=200, one sees that the median orientation size is 

if n < 2; 

1 if 2 < n < 20; 

2 if 20 < n < 42; 

3 if 42 < n < 67; 

4 if 67 < n < 93; 

5 if 93 < n < 120; 

6 if 120 < n < 147; 

7 if 147 < n < 174; 
,8 if 174 < n < 200. 

Note that by Theorem 2.6 there are exactly 8 orientations of size zero, 40 of size at most 1, 
104 of size at most 2, 232 of size at most 3, 488 of size at most 4, 1000 of size at most 5, 
2024 of size at most 6, 4072 of size at most 7 and 8168 of size at most 8. Thus for n = 27, 
just 104 orientations account for the majority of the tiles, while for n = 120, a mere 1000 
orientations account for the majority of the tiles. 

Table 1 shows the distribution of sizes for three samples, a microscopic sample with 
8 5 = 32768 tiles, a macroscopic sample with 8 27 ~ 2.42 x 10 24 tiles (about 4 moles), and an 
impossibly large sample of 8 120 ~ 2.35 x 10 108 tiles. Table 2 shows similar results for n = 10, 
20, 30, 40 and 50. 



4 Approach to the ergodic limit 

Let f(M) be a function on 5*0(3) and let f n (M) be the average of / over the 8 n tiles in the 
n-fold subdivision of a tile of orientation M. Equivalently, 



fn(M) = 8-» £ £ • • • £ /(M0 ilff 



■'2 



din)' 



As n — > oo, f n (g) approaches the constant limit J so ^ fdfj,, where d\i is Haar measure. The 
question is how fast the (ergodic) limit is approached. 

One approach to the problem is to consider the linear operator C : L 2 (SO (3)) — > 
L 2 (SO(3)) defined by 



Cf(M) 



1 



(4.1) 



i=i 



We can expand an arbitrary function in eigenfunctions ^ of C, with eigenvalues Xf 

f = E«^, ( 4 -2) 
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Size 


n = 8 


n = 27 


n = 120 





0.2588 


0.0971 


0.0145 


1 


0.5531 


0.3023 


0.0524 


2 


0.1719 


0.2931 


0.0788 


3 


0.0159 


0.1903 


0.1038 


4 


0.0003 


0.0847 


0.1221 


5 





0.0260 


0.1299 


6 





0.0055 


0.1256 


7 





0.0008 


0.1110 


8 





0.0001 


0.0897 


9 








0.0665 


10 








0.0453 


11 








0.0283 


12 








0.0163 


13 








0.0086 


14 








0.0042 


15 








0.0019 


16 








0.0008 


17 








0.0003 


18 








0.0001 


19 








0.0000 



Table 1: Distributions of Sizes in Various Samples. 



Size 


n= 10 


n = 20 


n = 30 


n = 40 


n = 50 





0.2199 


0.1272 


0.0878 


0.0654 


0.0508 


1 


0.5230 


0.3738 


0.2780 


0.2154 


0.1718 


2 


0.2199 


0.3023 


0.2842 


0.2488 


0.2140 


3 


0.0353 


0.1457 


0.2013 


0.2149 


0.2089 


4 


0.0019 


0.0426 


0.1011 


0.1422 


0.1635 


5 


0.0000 


0.0075 


0.0363 


0.0728 


0.1038 


6 





0.0008 


0.0093 


0.0290 


0.0537 


7 





0.0000 


0.0017 


0.0090 


0.0227 


8 





0.0000 


0.0002 


0.0021 


0.0079 


9 





0.0000 


0.0000 


0.0004 


0.0022 


10 





0.0000 


0.0000 


0.0001 


0.0005 



Table 2: Additional Distributions of Sizes. 
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so that 



/„ = J2 a i X i& ( 4 - 3 ) 

i 

Of course, the constant function is an eigenfunction of £ with eigenvalue 1. The rate at 
which /„ approaches a constant function is determined by the second-largest eigenvalue of 
£. 



The operator £ is quite similar to the Hecke operators considered in [|LPS1| , |LPS2 



except that the set {gi} is not balanced, as the inverses of and g$ are not in the set. As a 
result, £ is not a self-adjoint operator on L 2 '(SO '(3)), not even with a weighted measure, and 
we know of no a priori reason why the eigenvalues of £ need be real. However, the numerical 
experiments described below indicate that the eigenvalues of £ are indeed real. We do not 
understand the reason for this. 

The natural (left and right) actions of £0(3) on L 2 (SO(3)), namely (I, r) : f(g) — > 
f(lgr), decomposes L 2 (SO(3)) into irreducible representations of £0(3). Since £ is a lin- 
ear combination of right-translations, it acts separately on each representation. Thus the 
infinite-dimensional job of diagonalizing £ on all of L 2 breaks down into a sequence of 
finite-dimensional problems, namely diagonalizing £ on each irreducible representation that 
appears in the decomposition of L 2 . 

The irreducible representations of £0(3) are standard. The spin-£ representation, of 
dimension 2£ + 1, appears 21 + 1 times in the decomposition of L 2 . In particular, each 
representation appears at least once. For each representation Rg, we must diagonalize the 
(2£ + 1) x (2£ + 1) matrix 

C t = l[Re(gi) + --- + R e (gs)}. (4.4) 



We have done this, using MATLAB, for I up to 300, thereby diagonalizing £ on a 
subspace of L 2 (£0(3)) of total dimension greater than 36,000,000. There are numerous 
eigenvalues over 0.993, but none over 0.994; the largest known eigenvalue is A ~ 0.99381, 
appearing in the I = 258 representation. This proves Theorem 1.4, and suggests there 
is an upper bound to the eigenvalues slightly above 0.99381. Since powers of 2 appear 
repeatedly in the analysis, and since Theorem 1.2 suggests a time scale of 28 subdivisions, 
2~ 1 /( 4x28 ) ~ 0.9938303 is a fairly natural guess for that upper bound. 

Table 3 lists the largest eigenvalue for each of the first 20 representations. Table 4 lists 
the "record" eigenvalues, those larger than any eigenvalue corresponding to smaller £. The 
MATLAB program for finding the largest eigenvalues from the first Lmax representations is 

for L=l:Lmax 

% First specify angular momentum matrices 
jl=zeros(L+L+l) ; j2=jl; j3=j2; 
for m=l : (L+L) 

jl(m,m+l)=(L+L+l-m)*i/2; jl(m+l,m)=m*i/2; 

j2(m,m+l)=(L+L+l-m)/2; j2(m+l ,m)=-m/2; j3(m,m)=(L+l-m) *i ; 
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t 


Eigenvalue 


t 


Eigenvalue 


1 


0.50000 


11 


0.90142 


2 


0.78785 


12 


0.96793 


3 


0.50000 


13 


0.94705 


4 


0.92693 


14 


0.98797 


5 


0.62500 


15 


0.90546 


6 


0.94912 


16 


0.96437 


7 


0.84548 


17 


0.94949 


8 


0.98454 


18 


0.99048 


9 


0.73649 


19 


0.93382 


10 


0.95848 


20 


0.98042 



Table 3: Largest Eigenvalues for First 20 Representations 



i 


Eigenvalue 


1 


0.50000 


2 


0.78785 


4 


0.92693 


6 


0.94912 


8 


0.98454 


14 


0.98797 


18 


0.99048 


32 


0.99243 


45 


0.99324 


56 


0.99335 


72 


0.99362 


248 


0.99367 


258 


0.99381 



Table 4: Record Eigenvalues 
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end; 

j3(L+L+l,L+L+l)=-L*i; 

S = expm(pi*j2/2) ; T=expm(pi*j 1/3) ; % Exponentiate to get rotations 
A=S*T~3; B=S; C=S~2*T~3; D=S~2*T~2; E=T~4; 

F=(3*eye(L+L+l)+A+B+C+D+E)/8; G=eig(F) ; H(L,1)=L; H(L,2)=max(G) ; 
end; H 

5 Conclusions and Speculations 

Much of the interest in tilings stems from the fact that crystals and quasicrystals are modeled 
by tilings in which the tiles appear in only a finite number of orientations. Is there a form of 
matter that is modeled on aperiodic tilings with statistical rotational symmetry, in particular 
on a substitution tiling such as the quaquaversal? If so, how would such a form of matter 
be recognized? 

This paper provides the tools needed to recognize such a substance. We have computed 
the distribution of orientations in a sample of size 8™ of the QQ tiling. The exact same 
calculation gives the distribution of orientations of 8-tile clusters in a sample of size 8 n+1 , of 
64-tile clusters in a sample of size 8 n+2 , and so on. In particular, the angular distribution of 
the relative positions of nearby tiles, and therefore the angular distribution of the electron 
diffraction pattern, are governed by the results of this paper. 

To detect an asymptotically round quasicrystal, one must look at the diffraction pattern 
of a finite sample and decompose it into spherical harmonics. Some harmonics will be absent; 
those correspond to eigenfunctions of £ with eigenvalue A such that X n <^ 1. Other modes, 
with 1 — A comparable to or less than 1/n, will appear, however. A macroscopic (say, 
N = 8 27 ) sample of the QQ tiling, for example, would have a diffraction pattern with a small 
but probably detectable £ = 4 component, and with essentially undamped £ = 6, £ = 8 and 
£ = 10 components. There are relatively few eigenfunctions with such large eigenvalues (even 
for £ = 8, all but two of the 17 eigenvalues are less than 0.65), so the diffraction pattern 
will be much more isotropic than that of a crystal or ordinary quasicrystal, while much less 
isotropic than that of an amorphous substance such as a glass. 

The radial dependence of the diffraction pattern also carries important information. As 
with other quasicrystals, and unlike amorphous media, one expects asymptotically round 
quasicrystals to have radially self-similar diffraction patterns. In combination, a self-similar 
radial pattern combined with a mildly anisotropic angular pattern should indicate an asymp- 
totically round quasicrystal. 

We thank Charles Radin for helpful discussions and Margaret Combs for her Tj^Xnical 
expertise. B.D. thanks J.M. Linhart for her guidance and support. The work of L.S. is 
partially supported by NSF Grant No. DMS-9626698 and Texas ARP Grant 003658-152. 



16 



References 



[CR] J. Conway and C. Radin. Quaquaversal tilings and rotations, Inventiones Math. 132 
(1998), 179-188. 

[LPS1] A. Lubotsky, R. Phillips and P. Sarnak. Hecke Operators and Distributing Points 
on the Sphere I. Comm. Pure Appl. Math. 39 (1986) S149-S186. 

[LPS2] A. Lubotsky, R. Phillips and P. Sarnak. Hecke Operators and Distributing Points 
on 5 2 II. Comm. Pure Appl. Math. 40 (1987) 401-420. 

[Rl] C. Radin. Symmetry and Tilings, Notices Amer. Math. Soc. 42 (1995), 26-31. 

[R2] C. Radin, private communication. 

[RSI] C. Radin and L. Sadun. Subgroups of 5*0(3) associated wtih tilings, J. Algebra 202 
(1998), 611-633 

[RS2] C. Radin and L. Sadun. On 2-generator subgroups of 5*0(3), Trans. Amer. Math. 
Soc, to appear. 

[S] F. Su, Convergence of random walks on the circle generated by an irrational rotation. 
Trans. Amer. Math. Soc. 350 (1998), 3717-3741. 



17 



